Reference Files/Mann_Whitney_Stats.R

#' Mann-Whitney U Stats
#
#' This function takes two (2) data vectors as input
#' and outputs U statistics

#
# Some useful keyboard shortcuts for package authoring:
#
#   Install Package:           'Ctrl + Shift + B'
#   Check Package:             'Ctrl + Shift + E'
#   Test Package:              'Ctrl + Shift + T'

#' @param A A vector of values for variable A
#' @param B A vector of values for variable B

#' @return UA Value of UA
#' @return UB Value of UB
#' @return U Mann-Whitney U statistic (the lesser of UA and UB)
#' @references Chechile, R.A. (2020). Bayesian Statistics for Experimental Scientists. Cambridge: MIT Press

Mann_Whitney_U_Stats<-function(A, B){
  UA_vector<-rep(NA, length(A)) # UA counter
  UB_vector<-rep(NA, length(B)) # UB counter
  for (i in 1:length(A)){
    UA_vector[i]<-sum(A[i]>B)
  }
  for (j in 1:length(B)){
    UB_vector[j]<-sum(B[j]>A)
  }

  list(UA=sum(UA_vector),
       UB=sum(UB_vector),
       U = min(UA, UB))
}
danbarch/dfba documentation built on Jan. 30, 2024, 6:51 p.m.